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Abstract 

Here we report the polarization dependent non-resonant second harmonic generation (SHG) mea- 
surement of the interfacial water molecules at the aqueous solution of the following salts: NaF, 
NaCl, NaBr, KF, KCl, and KBr. Through quantitative polarization analysis of the SHG data, 
the orientational parameter 'D,(D={cos6)/{cos^6)) value and the relative surface density of the 
interfacial water molecules at these aqueous solution surfaces were determined. From these results 
we found that addition of each of the six salts caused increase of the thickness of the interfacial 
water layer at the surfaces to a certain extent. Noticeably, both the cations and the anions con- 
tributed to the changes, and the abilities to increase the thickness of the interfacial water layer 
were in the following order: KBr > NaBr > KCl > NaCl ^ NaF > KF. Since these changes can 
not be factorized into individual anion and cation contributions, there are possible ion pairing or 
association effects, especially for the NaF case. We also found that the orientational parameter D 
values of the interfacial water molecules changed to opposite directions for the aqueous solutions 
of the three sodium salts versus the aqueous solutions of the three potassium salts. These findings 
clearly indicated unexpected specific Na"*" and K"*" cation efi^ects at the aqueous solution surface. 
These effects were not anticipated from the recent molecular dynamics (MD) simulation results, 
which concluded that the Na+ and K+ cations can be treated as small non-polarizable hard ions 
and they arc repelled from the aqueous interfaces. These results suggest that the electrolyte aque- 
ous solution surfaces are more complex than the currently prevalent theoretical and experimental 
understandings. 
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I. INTRODUCTION 



Recently, aqueous solution surfaces of simple electrolytes have attracted much attention 
for its importance in the broad chemical, environmental and biological processes at the 
molecular interfaces and membranes. Conventionally, it has been widely 

believed that the aqueous solution surfaces were empty of ions. This classical picture was 
based on the macroscopic surface tension measurement of the electrolyte solutions, and was 
theoretically described by the negative excess of the ions from the thermodynamic Gibbs 
adsorption equation, as well as the microscopic image charge repulsion theory. 7, ^ Recently, 
new findings in the atmospheric chemical research, direct surface measurement techniques 
and molecular dynamics (MD) simulations have challenged this classical picture of ions 



depletion at the electrolyte solution surface region since Gibbs and Onsager. 
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Consequently, theoretical simulations and experimental techniques have been employed 
to provide answers to the two equally important questions about the electrolyte aqueous 
solution surface. One is how the anions and cations are spatially distributed in the surface 
region, and the other is how the interfacial water molecules are affected by this ion distri- 
bution. To answer the former question requires techniques which can directly measure the 
presence and profile the depth of the ion distribution at the liquid surface; while to answer 
the later requires techniques which can directly measure the presence and change of the in- 
terfacial water molecules. The convergence of and the correlation between the two answers 
may provide a through picture for the understanding of the electrolyte aqueous solution 
surface. The fact is that the conclusions of most of the recent theoretical and experimental 
studies converged on the picture that the anions, especially the larger and more polarizable 
anions, are presented at the interfacial region of various simple inorganic salt solutions. 

Molecular dynamics (MD) simulation results using the polarizable force fields concluded 



the anions, such as the larger and more polarizable I and Br 
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18( 1 MD simulation also revealed that the anions and 



for significant enrichment^o. 
at the interfacial region, 
cations should be segregated for the electrolytes consist of more polarizable anions at the 
electrolyte aqueous solution surface, with the double layer structure in which the anions 
enriched at the top surface layer and the cations resided below. [5] Except for some voices of 



discontent. 
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20| the majority of the experimental studies on the electrolyte aqueous so- 
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lution surfaces, including the studies with X-ray photoelectron spectroscopy (XPS), 
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X-ray fluorescence 



211 ] and nonlinear optical 



spectroscopy 



second harmonic generation (resonant SHG) 
eration vibrational spectroscopy, (SFG-VS) 
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techniques, such as the resonant 
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3G, 
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2m and sum-frequency gen- 
32l | lent their support to this 



new physical picture with the more polarizable anions enriched at the electrolyte aqueous 
solution surface. 

Since in this new picture the more polarizable anions are expected to be present at the 
electrolyte aqueous solution surface, while the small non-polarizable hard cations, such as 
Na^ and K"^ etc., and anions, such as F~, are expected to be repelled from the aqueous 
solution surface, experimental studies have been focused mainly on either to prove the above 



conclusions 
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on the interfacial water molecules. 
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33l or to investigate the specific anion effects 
271 . 134| Therefore, to our best knowledge, so 



far there has been no systematic investigation on the cation effects on the interfacial water 
molecules at the electrolyte aqueous solution surfaces. [3|, l5|, IgJ 

Recently, we developed and employed the surface sensitive non-resonant SHG technique 



to measure the polarization dependent SHG response 



Tom the interfacial water molecules 



at the NaF, NaCl, and NaBr salt solution surfaces. [141. |35[| In this study, in order to ensure 
measurement of small changes in the SHG signal, we first developed effective procedures 
to monitor and remove the impurities in the salt solution samples. Then, the quantitative 
polarization analysis of the measured SHG data showed that the average orientation of the 
interfacial water molecules changed slightly with the increase of the bulk concentration of the 
NaF, NaCl and NaBr salts from that of the neat air/water interface. The observed significant 
SHG signal increase with the bulk salt concentration was attributed to the overall increase 
of the thickness of the interfacial water molecular layer, following the order of NaBr > NaCl 
~ NaF. The absence of the electric-field-induced SHG (EFISHG) effect indicated that the 
electric double layer at the salt aqueous solution surface is much weaker than that predicted 
from the molecular dynamics (MD) simulations. These results provided quantitative data 
to the specific anion effects on the interfacial water molecules of the electrolyte aqueous 
solution, not only for the larger and more polarizable Br~ anion, but also for the smaller 
and less polarizable F~ and Cl~ anions. 

The most intriguing thing in this study was that we were able to quantitatively measure 
weak concentration dependent F~ anion effect of the NaF salt on the interfacial water 
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molecules, {l^ even thoiigh this F~ anion effect was expected to be insignificant from the 
recent MD simulations. |3|, l33| This immediately led us to look at whether the F~ anion 
effect of the KF salt aqueous solution surface can also be observed. Moreover, there is also 
one known advantage to investigate the KF aqueous solution because this can extend the 
bulk concentration range of the F~ anion. The saturation solubility for NaF in water at 
20°C is only 0.98M (4.13g NaF per lOOg water), while the saturation concentration in KF is 
15. 4M (89. 8g KF per lOOg water). 36j Surprisingly, we not only observed the significant F~ 
anion effect of the KF salt aqueous solution surface up to much higher bulk concentration as 
expected, we also found that the orientational parameter of the interfacial water molecules 
went to the opposite direction to that of the NaF aqueous solution surface as the bulk 
concentration of NaF and KF salt increase. Therefore, we started doing experiment also on 
the surface of the KCl and KBr salt aqueous solutions, and we observed distinctive specific 
Na^ and K+ cation effects on the interfacial water molecules at the NaF, NaCl, NaBr, KF, 
KCl, and KBr salt solution surfaces. 

In the following sections, we shall report the non-resonant SHG measurement and analysis 
of the specific Na"*" and K+ cation effects on the interfacial water molecules at these six 
solution surfaces. We shall show that while the thickness of the interfacial water layer at 
the surfaces of all the six salt solutions increases as the bulk concentration increase, the 
orientational parameter D values of the interfacial water molecules changed to opposite 
directions for the aqueous solutions of the three sodium salts versus the aqueous solutions of 
the three potassium salts. Moreover, the abilities to increase the thickness of the interfacial 
water layer were found to be in the following order: KBr > NaBr > KCl > NaCl ~ NaF > 
KF. These specific Na"*" and K"^ cation effects, as well as the F~, Cl~ and Br~ anion effects, 
indicate that the electrolyte aqueous solution surfaces are likely to be more complex than 
what the currently prevalent theoretical and experimental understandings can offer. 

II. THEORETICAL BACKGROUND FOR THE SHG DATA TREATMENT 

The theoretical tools used in this report shall fully follow the detailed description in the 



previous published paper, since this work is an extension from the previous study. IJ] Here 



we only present what is necessary to understand the analysis of the SHG data in this report 



The SHG intensity I{2uj) reflected from the electrolyte aqueous interface is given by 37|, 



5 



38|,l39| 

327r^ti;^ sec^ /? , ,2 n / n 

Here /(cu) is the incoming laser intensity, Cq is the speed of the hght in the vacuum, (3 is 
the incident angle from the surface normal, defined as z axis in the laboratory coordinates 
system (x, ni{uji) is the refractive index at the frequency uoi of the medium in which 
the laser beam propagates. Xe// is the effective macroscopic second order susceptibility. 
In the SHG experiment, because there are only one incident laser beam and one out-going 
signal beam, there are three independent polarization combinations for the Xeff'- namely the 
s-in/p-out {xeff,sp), the 45°-in/s-out (Xe//,45°s) and the p-in/p-out {Xeff,pp)- Here p denotes 
the polarization in the incident plane, while s the polarization perpendicular to the incident 



plane. With the microscopic loca. 
Fresnel factor Las, the Xe//S are 
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d factors incorporated implicitly into the tensorial 



Xeff,sp = L,,{2uj)Lly{uj)sin/3xzyy 
XeffA5°s = Lyy{2u)L,,{uj)Lyy{u)sinpXyzy 

Xeff,pp = +L,z{2uj)Ll^{uj)sinf3cos'^f3xzxx 

-2Lx^{2u)L^^{uj)L^^{u)sinPcos^PXxzx 
+L,,{2uj)Ll^{uj)sin^(3xzzz (2.2) 

in which the XijkS are the seven non-zero susceptibility tensors, i.e. Xzzz, Xzxx = Xzyy, 
Xxzx = Xxxz = Xyzy = Xyyz, of the rotatioually symmetric interface along the interface 
normal z, with x in the incident plane. 

In the SHG measurement, the polarization is usually fixed at either p ot s, while the 
polarization angle G of the incident laser beam is varied for the full 360° using a half 



waveplate. Then, one has 35|, |40| 



h ^ IXe//,pP = |Xe//,ppCOS^0 + Xe//,spSin2 0P 

Is oc |xe//,sP = |Xe//,45°sSin29p (2.3) 

In the case of non-resonant SHG, each of the Xeff,pp, Xeff,sp and Xe//,45°s terms is real, 
and their relative values can be obtained from fitting the experimental data with the Eq J2.3l 



by knowing the relative signs between them. The relative signs can be determined with 



the measurement of the SH signal in the 45° polarization as described elsewhere. [351] It was 
also shown that using the Eqs l2.3[ Xeff,pp, Xeff,sp and XefiA5°s terms can be most accurately 



35| 



determined from the polarization dependent SHG data. 

In the previous study, we showed that the contributions from the ions and from the electric 
field induced SHG were negligible for the non resonant SHG signal from the electrolyte aque- 
ous solution surface. Therefore, only the contribution from the interfacial water molecules 
need to be considered. For molecules with the symmetry, such as the water molecule, 
there are seven nonzero microscopic polarizability tensor elements f]ijk{i,j,k = a,b,c) in 
molecular coordinates system (a,6, c): namely, Pccc, Pcaa, Pcbb, Paca, Paac, Pbcb and Pbbc, with 
Paca = Paac and Pbcb = Pbbc fof the SHG process. When the optical frequency is far below 
the resonance, (3ccc is usually negligible, by defining R = {(3caa + Pcbb)/{Paca + Pbcb), one 
has,Q,Q 



Xzx. {R-2)D + {R + 2) 

(2.4) 



X,z. -RD + (i? + 2) 

Xzzz _ 2{R + 2)D-2{R + 2) 



X.zx -RD + iR + 2) 

Here D = (cosd) / {cos^O) is the orientational parameter with 9 as the tilt angle of the 
water dipole axis from the interface normal, () denotes the ensemble average over the whole 
orientational distribution. Because the two ratios XzxxIXxzx and XzzzlXxzx can be accurately 
obtained from the polarization SHG measurement^the R and D values can be readily 



obtained with good accuracy, using the Eqs J2.4[ |35l . |40 1 

Now with the orientational parameter D known, the relative number of water molecules 
contributing to the measured SHG signal from the interface of the different electrolytes 
as well as different electrolyte concentrations can be calculated from the changes of the 
corresponding SHG intensities. Simply to put it, 

hu. oc \xeff? oc \N,Xefffl? = A^'Ir^T (2.5) 

Here A^^ is the interfacial number density of the water molecule, Xeff,o is the average 
effective susceptibility per molecule. Here Xe//,o can also be written into the per molecule 
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orientational functional r{6), which is directly related to the orientational parameter D. 35|, 



III. EXPERIMENTAL DETAILS 



The procedures to pretreat the salt sample and purify the salt solution were described 
in detail previously and should be referred when in doubt. 1^ The necessary experimental 
detail and the information for sample preparation are described below. 



A. SHG measurement 



The SHG experimental setup is as in the previous reports. [1^, [371, |40(] A broadband 
tunable mode-locked femtosecond Ti:Sapphire laser (Tsunami 3960C, Spectra-Physics) is 
used for the reflected-geometry SHG measurement. Its high-repetition rate (82 MHz) and 
short pulse width (80 fs) make it suitable for detection of the weak second-harmonic signals. 
Its long term power and pulse-width stability also make it easy for quantitative analysis of 



the SHG data. 4l| The 800 nm fundamental laser beam is focused on the solution interface 
at the incident angle of (3=70°, and the SHG signal at 400 nm is detected with a high-gain 
photomultiplier tube for single photon counting (R585, Hamamatsu) and a photon counter 
(SR400, Stanford Research Systems). Typically the dark noise level is less than 1 counts/s, 
bett. than .ported p.ev,ous.,Q TKe typicanase. power is 500 mW. The efficiency of 
the detection system for the p polarization is 1.23 times of that for the s polarization. A 
cylindrical Teflon beaker (050x6mm) is used to host the solution. The optical polarization 
control and the SHG data acquisition are programmed and collected with a PC. It takes 
about 150s to 300s to collect the polarization curve. The room temperature is controlled at 
22.0±1.0°C, and the humidity in the room is controlled around 40%. 



B. Sample preparation 

In the experiments, the liquid water was purified with a Millipore Simplicity 185 
(18.2Mf2-cm) from double distilled water. The KF, KCl and KBr used were purchased 
from the Sigma- Aldrich (ACS reagent grade, 99%+ purity). Before the preparation of the 
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salt solution, the salts were baked at around 500°C for more than 6 hours in order to remove 
the organic impurities. The glassware was cleaned with hot chromic acid and then rinsed 
with Millipore water for several times. The solutions were filtered with a syringe with a 
0.22yum membrane filter made of PVDF material by Millipore (Durapore GVWP01300) to 
remove the possible insoluble particles existing in the solution. The SFG-VS spectral and 
SHG fluctuation measurement are used to monitor the contaminates. 



IV. RESULTS AND DISCUSSION 

A. Quantitative analysis of the polarization dependent non-resonant SHG data 

In the previous study, jl^ we showed that the non-resonant SHG signal from the aqueous 
solution surface of simple electrolytes was originated from the surface water molecules, and 
the contribution from the surface ions was negligible. We also showed that the contribution 
to the non-resonant SHG signal through the electric field induced second harmonic gener- 
ation (EFISHG) mechanism possibly from the surface potential due to the surface charge 
separation was negligible. Careful sample treatment and purification procedures were also 
developed to ensure data reproducibility. Here we just followed the experimental and data 
analysis procedures established in the previous literature, from which the accuracy and re- 
producibility of the non-resonant SHG data and the quantitative analysis of the SHG data 
were assured. Q, 

Fig. [1] shows the non-resonant SHG intensities for the NaF, KF, NaCl, KCl, NaBr and 
KBr solution surfaces in the three independent polarization combinations, i.e. p-p, 45°-s 
and s-p polarization combinations, at different concentrations. The data for the NaF, NaCl 



and NaBr were published in our previous work, [ij] The saturation concentrations of KF, 
KCl and KBr at the 20°C are 15.4M (89.8g KF per lOOg water), 4.6M (34.2 g KCl per lOOg 
water), and 6.3M (65. 3g KBr per lOOg water), respectively. ^36] In comparison, the saturation 
concentration of the NaF, NaCl and NaBr at 20°C are 0.98M(4.13 g NaF per 100 g water), 
6.2 M (36.0 g NaCl per 100 g water), and 9.2 M (94.6 g NaBr per 100 g water). [36;] As 
mentioned above, the large saturation concentration for KF enabled us to go to higher bulk 
concentration with the F~ anion by using the KF salt; while with the NaF, the saturation 
concentration is only 0.98M. In the experiment, the highest bulk concentration for KF, KCl, 



9 



and KBr are 5M, 4M and 4M, respectively. 

It is clear from the Figure [1] that for all the salts, the SHG signal increased almost 
linearly with the increase of the bulk electrolyte concentration, and the general behaviors 
for the potassium salts are similar to those of the sodium salts, with the biggest slope for 
the KBr and smallest slope for the KF. Therefore, this suggests that similar to the NaF, 
NaCl and NaBr salts, Q, KF, KCl and KBr salts also resulted in the increase of the 
thickness of the interfacial water molecular layer. However, a closer examination of the SHG 
data revealed apparent differences between the NaF/KF, NaCl/KCl and NaBr/KBr pairs, 
indicating non-negligible specific Na"*" and K"^ cation effects on the non-resonant SHG signal 
of the interfacial water molecules. For example, in the Figure [T], the slopes of the increase 
of the p-p signal for the KF, KCl and KBr are larger than the slope for the NaF, NaCl and 
NaBr, respectively; while for the 45°-s and s-p data showed completely the opposite trend. 
This relative differences between the intensities of different polarization combinations usually 
indicate possible differences of orientational order of the surface molecules contributing to 



the SHG signal.|35|, 133,1381 



The differences of the orientational order parameter D values are clearly presented in the 
Figure [21 The D values were obtained by using Eqs. 12.21 and Eqs. 12. 4[ and the calculated 
results for the KF, KCl and KBr solutions are listed in the table HI In the table [U the cal- 
culated molecular polarizability ratio R = {Pcaa + Pcbb) / {Paca + Pbcb) values for all three salts 
at different bulk concentrations were all the same within the small error bar as the value 



calculated from the SHG data for the neat air/water interface. 40|] These values are also the 



same as those obtained for the NaF, NaCl and NaBr aqueous solution surfaces as reported 
previously, [l^ This further confirmed that the assumption that the non-resonant SHG ob- 
served from the KF, KCl and KBr aqueous solution surfaces were indeed from the interfacial 
water molecules. The D and 6 values for the NaF, NaCl and NaBr aqueous solution surfaces, 

n 

which were published previously. |14I] are also plotted in the Figure [2] for comparison. Even 
though the D values only changed slightly at different bulk salt concentration, the difference 
is clearly bigger than the error bar as the bulk salt concentration increase. It is interesting 
to see that the D values for the potassium salt series are apparently larger than the D values 
for the sodium salt series. This also confirms the qualitative surmise above based on the 
different trends for the polarization dependence of the non-resonant SHG signal for the KX 
(X=F, CI or Br) salt series and the NaX salt series. 
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In the Figure [21 the 6 value is the apparent average tilt angle between the dipole vector of 
water molecule and the surface normal calculated using the assumption of a 5 distribution 
function for this tilt angle. The assumption for a 5 distribution function of 9 is certainly not 
correct. However, as discussed intensively in the literature, because within a certain range 
of the width for the orientational distribution, this 9 value can be used to represent the 
center value for the orientational distribution and can be used to simulate the orientational 
dependent total SHG susceptibility, i.e. the orientational functional r(^).[37|, |42|, |43| This 
point shall be further illustrated below. 

With the orientational parameter known, the relative surface density Ng of the interfacial 
water molecules which contribute to the measured non-resonant SHG signal can be readily 
calculated using the Eqs. 12.51 These values for the KX salts are also listed in the Table HI 
Figure [3] shows the Ns values for the NaF, KF, NaCl, KCl, NaBr, and KBr salt aqueous 
solution surfaces. The Ng values for the three NaX salts are the same as from the previously 
published results, [l^ It is clear the Ng for each salt increases almost linearly with the bulk 
salt concentration. As discussed previously, l4] since the probing area of the laser is the same, 
the increase of the Ng represents the increase of the average depth of the water molecules 
in the interfacial region which contributes to the non-resonant SHG signal. Now we can 
fit the Ng versus the bulk concentration of each of the salts data with a linear function of 
Ng = Ngfi + kC. Here A^^.o = 1 denotes the surface density of the neat air/water interface, C 
is the bulk salt concentration, and the slope k represents the ability to increase the thickness 
of the surface water molecules. The fitting results of k values are listed in the Table [Tll The 
unit for k is M~^. 

According to the k values in the Table [Tll the abilities to increase the surface water 
thickness for the six salts follow the following order: KBr > NaBr > KCl > NaCl ~ NaF > 
KF. Clearly, both the anions and cations contributes to the increase of A^^- In the previous 
study with the NaF, NaCl and NaBr salts, the order was NaBr > NaCl ~ NaF, and mostly 
the effect was attributed to the ability to be polarized for the larger anions. However, here 
we clearly have KBr > KCl > KF, and one can also conclude that even though there are 
differences on the abilities to increase the aqueous surface thickness for the individual ions, 
the abilities for the cations and anions are generally with comparable order of magnitudes. 
Therefore, the cation effects can no longer be ignored. These findings on the specific anion 
and cation effects from the quantitative analysis of the non-resonant SHG data imply that 
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the electrolyte aqueous solution surface can be more complex than previously thought. 



B. The structure and reorientation of the interfacial water molecules at the elec- 
trolyte solution surfaces 

There have been many theoretical and experimental studies on the structure of the in- 
terfacial water molecules. The vibrational spectral feature around 3700cm~^ of the free 
0-H stretching vibration was first observed by Shen and co-workers in 1993 using the sur- 
face specific sum frequency generation vibrational spectroscopy (SFG-VS).|44| Subsequent 
SFG-VS experimental studies and quantitative analysis further concluded that this water 
molecule with the free 0-H straddles across the air/water interface with its dipole vector 
almost parallel to the interface, i.e. on average the free 0-H tilted from the surface normal 



around 30-35°. 
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471, 
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491] This latter conclusion was also su ppo rted b y molecular 



51 



52 



53, 
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571, 



58, 



59 



60| 



dynamics (MD) and Monte Carlo (MC) simulations. [50 
SFG-VS experiment also revealed that the SFG spectra in the 3200-3600cm~^ region were 
from the contributions of the hydrogen bonded water molecules slightly below the topmost 
water layer at the air/water interface, and these water molecules had both of their hydro- 
gen atoms hydrogen bonded. 4^, 47] These are the non-straddle-type water molecules at the 
air/water interface, and the SHG signal from the air/water interface comes mainly from the 
non-straddle-type water molecules. [4o| 

The SFG-VS studies by Allen and Richmond and their co-workers on the NaF, NaCl, 
NaBr and Nal aqueous solution surfaces concluded that the free 0-H spectra was almost 
undisturbed with the addition of the salts, while the hydrogen bonded water spectra were 
variously affected. 19|, l27|] Using the more informative phase- sensitive sum-frequency spec- 
troscopy technique, Shen et al. recently investigated the hydrated I~ anion effects on the 
interfacial water structure. 3J] Their conclusion was that the presence of I~ anion near the 



interface region may not disturb the water molecular structure at the topmost surface layer, 
while it can reorient the water molecular in the subphase. These conclusions are consistent 
with our non-resonant SHG results on the interfacial water molecules at these electrolyte 
aqueous solution surfaces. Moreover, our SHG results provides not only quantitative mea- 
surement on the slight reorientation of the interfacial water molecules, but also quantitative 
measurement on the increase of the average thickness of the interfacial water layer, for the 
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electrolyte aqueous solution surfaces. 

Here we would like to further discuss the different reorientation effects on the interfacial 
water molecules by the different anions and cations. 

Firstly, the change of the values of the orientational parameter D in the whole bulk 
concentration range is small. If we do not compare the results for the KX series with the 
NaX series, we may conclude that there is essentially no apparent change of the D values 
for the KX series. This indicates that the reorientation effect for the KX series is rather 
insignificant. Therefore, the significant increase of the SHG signal is essentially not the 
result of water molecule reorientation. 

Secondly, from the Figure [21 the direction of the reorientation effect is mainly determined 
by the Na"*" and K+ cations, instead of the anions. This finding is intriguing because neither 
the classical picture nor the recently revised picture can be used to predict this significant 
cation effect. In the classical theory, both the anion and cations are thought to be repelled 
from the aqueous solution surface. The recent MD simulations which predicted enrichment 
of more polarizable anions also negated the presence of the non-polarizable cations, as well 
as the non-polarizable anions, at the aqueous solution surface, The results here suggested 
that probably even for the much less polarizable cations and anions can still be present or 
enriched at the interface region. It is probable that in order to keep the neutrality of the anion 
enriched interface region, more cations are attracted to the surface region. This picture is 
certainly consistent with the fact we reported previously that there was insignificant EFISHG 



contribution to the total SHG signal. 



14| However, in order to find reconciliation with the 



current picture of charge segregation from the MD simulation results, much effort are still 
needed. 

With the known R values, the orientation dependent per molecule Xe//,o values in different 
polarization combinations can be simulated using the Eqs J2.2l and Eqs J2.4[ When the Ns 
value was known, the Xe//,o value relative to the neat air/water interface value Xeff, water can 
be calculated with the Eqs J2.5[ Comparison of these two values can help us understand how 
the SHG signal changes with the averaged molecular orientation of the water molecules at 
the electrolyte aqueous solution surfaces. Figure H] shows the simulation results of the Xe//,o 
values in the p-p, 45°-s and s-p polarization combinations using R=0.69 and {Paca + Pbcb)=^- 
From the curves in the Figure HI we can see that if the average orientation angle 6 becomes 
slightly larger than the neat air/water interface value ^=40.1°, the Xpp,o value shall become 
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larger, while the X45°,o and Xsp,o values shall become smaller. This is indeed the case for 
the KF, KCl and KBr aqueous solution surfaces from the SHG measurement. On the other 
hand, if the average orientation angle 6 becomes slightly smaller than the neat air/water 
interface value 6'=40.1°, the Xpp,o value shall become smaller, while the X45°,o and Xspfi values 
shall become larger. This latter case is consistent with the SHG data for the NaF, NaCl 
and NaBr aqueous solution surfaces. 

In the Figure m we also presented the simulation of the Xe//,o values when the orienta- 
tional distribution function is not the 6 function. It is clear that when the full width at 
the half maximum (FWHM) of the orientational distribution function is 30°, the trend of 
dependent change of each Xeff,o component remained the same in the vicinity of 6'=40.1°. 
For example, just like the 6 distribution case, when the average orientation angle 6 be- 
comes slightly larger than the neat air/water interface value 6'=40.1°, the Xppfi value shall 
become larger, while the X45°,o and Xsp,o values shall become smaller. Our simulation further 
demonstrated that as long as FWHM does not exceed 65°, this trend shall remain the same. 
Therefore, the conclusions that the NaX salt series slightly reorient the interfacial water 
molecules to smaller average tilt angle and the KX salt series slightly reorient the interfacial 
water molecules to larger average tilt angle remained unaffected when the orientation of 
the interfacial water molecules possess a broad range of distribution width. This is another 
example to illustrate that in the polarization and orientation analysis of SHG, as well as 
SFG, the orientational order can be understood even though the orientational distribution 
is in a broad range. 42f| Thus, the small shift of the average orientation of the interfacial 
molecules, as representatively shown in the Figure [5l can be accurately determined through 
polarization and orientational analysis of the SHG/SFG data. 

The unexpected specific cation effect on the reorientation of the interfacial water 
molecules by the Na^ and cations as determined here indicates complex interactions 
between the ions and the water molecules at the electrolyte aqueous solution surfaces. It is 
rather surprising to see that even though the F~, Cl~ and Br~ anions influenced the extent 
of the reorientation effect, they do not control the direction of the reorientation. These are 
surprising results from our current understandings on the ion adsorption to the aqueous 
surfaces. 

Based on these different cation effects on the hydrogen bonded water molecules in the 
surface region, we surmised that there might be some observable effects on the hydrogen 
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bonding 0-H vibrational spectra for these interfacial water molecules. This prompted us to 
measure the SFG vibrational spectra of the interfacial water molecules in the 3000-3800cm~^ 
region at the NaF and KF aqueous solution surfaces. The SFG-VS spectra of the NaF 
aqueous solution surfaces have been investigated by various research groups. However, these 
data were not quite consistent with each other. [2], 14, 19, 27| Moreover, to our knowledge. 



there has been no SFG-VS measurement in the literature on the water molecules at the 
KF aqueous solution surface, partly because the cation effects between the Na"*" and K+ 
on the surface water molecules were expected to be negligible in the best of the current 



understandings. 
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33| 



We measured and quantitatively analyzed the SFG-VS spectra of the neat air /water in 
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61| Based on these successful 



terface and other water surfaces in the past few years, 
experiences, we measured the SFG-VS spectra of the NaF and KF aqueous solution sur- 
faces at different bulk concentrations. We found that the NaF water SFG-VS spectra were 
consistent with the reported spectra by the Richmond group, Q, Q and we also observed 
distinctive differences in the hydrogen bonding 0-H stretching region in 3000-3600cm~^ for 
the NaF and KF aqueous solution surfaces. These results firmly confirmed the specific 
cation effects between the Na"*" and K"*" as measured with non-resonant SHG in this work. 
The specific ion effects on the orientation of the free 0-H group of the straddle type surface 
water can be more subtle. Since SHG is not sensitive to the straddle type of surface water 
molecules, further polarization dependent SFG-VS spectra measurement and analysis have 
been conducted. The SFG-VS results on the specific cation effects, especially the Na"*" and 
K"*" cations, shall be reported elsewhere. [62l] 



C. Complexities of the specific anion and cation effects on the interfacial water 
molecules 



With the data on the NaF, NaCl, NaBr, KF, KCl, and KBr solutions, and with both the 
specific anion and cation effects established, more detailed comparison between the anion 
and cation effects on the interfacial water molecules can be explored. 

In the bulk solution, the anion and the cation interactions with water molecules can 



influence the macroscopic properties of water. |63l. l64| such as the surface tension [651] and 
the viscosity. 66f] The hydrated ions effects on the hydrogen bonded networks in the bulk 
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water are the key questions people want to understand. [63| For a long time, the qualitative 
Hofmeister series first proposed more than a century ago were used to explain the ion 
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m It has 



effects on forming and breaking of hydrogen bonds in the aqueous phase, 
been believed that anions have more pronounced effects on the water structure than the 
cations. jsS] In the mean time, different cations can also influence the structure of water 
molecules. For example, the different behaviors of Na"*" and K"*" cations play important 
roles in many biochemical processes, such as in the ion channels crossing the biological 



membranes. Studies showed that these different 
hydration number and hydration free energies. [69 



aehaviors were largely due to their different 
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721] Recently, using SFG vibrational 



spectroscopy, Cremer and co-workers investigated the ion effects on the water molecules at 



the polymer surfactant aqueous interfaces. 



7J] In these studies, significant anion effects 



and insignificant cation effect were observed. 

From the Hofmeister series in the bulk, the Na"*" cation was traditionally called the 
'structure making' cation, while the K"*" cation was considered 'structure breaking' to some 
extent. 651] Using the mixed quantum mechanical (QM) and molecular mechanical (MM) 
simulation, Tongraar et al. reported that the Na"*" cation can order the structure of solvent 
molecules beyond the first hydration shell, while K"*" cannot go beyond the first hydration 
shell. [3] So it was argued that the Na"*" be recognized as the 'structure-making', while K"*" 
be 'structure-breaking'. Mancinelli et al. showed that water molecules are more orientation- 
a... disordered ..oud the ion than the Na^ ion.Q On the other hand. Saykahy . al. 
used the X-ray absorption spectroscopy and showed that the dissolved monovalent cations 
like Na"*" and K"*" did not have significantly different effects on the nature of hydrogen bonds 
around these ions. [77| Otma et al. found, using the femtosecond pump-probe spectroscopy. 



that the presence of ions does not lead to an enhancement or a breakdown of the hydrogen- 
bond network in liquid water, [t^ There is still no consensus regarding the ions effects on 
the hydrogen bond network of the water molecules. 

In the non-resonant SHG results, the KX series slightly increased the orientational pa- 
rameter D value of the interfacial water molecules, while the NaX series slightly decreased 
the D value. The larger its averaged D value, the more possibly tilted or disordered the in- 
terfacial water molecules. This seems to suggest that the Na"*" cation is 'structure-making', 
while the K"*" is 'structure-breaking'. However, in the bulk Hofmeister series, F~ was usu- 
ally considered 'structure making', while the Cl~ and Br~ 'structure breaking'. If judged 
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with the change of D values, this seemed to work for the KX salt series, but totally to the 
opposite for the NaX series, as shown in Figure [2l Therefore, we can only conclude that the 
ion effects at the electrolyte aqueous solution surfaces as observed in the non-resonant SHG 
measurement may be more complex than the Hofmeister series can qualitatively explain. 

One example of the possible complexities is that further examination of the k (slope) 
values in the Table [TTl suggested that the anions and cations may associate in some cases, 
especially for the NaF aqueous solution surface. If there is no anion and cation association at 
the interface, the k values in the Table [ITl should be able to factorize into ion specific additive 
factors. This can be achieved easily for the k values for the KF, KCl, KBr, NaBr and NaCl 
without considering the k value for the NaF. NaF stands out not only in this, but also in its 
saturation concentration (less than IM), which is far smaller than the other five salts (all 
larger than 4.5M or more). It has been long known that the F~ is significantly different from 
the larger Cl~ and Br~ ions in terms of interactions with the water molecule. js] However, 
there are huge differences between the saturation solubility of the LiF (O.IOM, or 0.27g LiF 
per lOOg water), NaF(0.98M, or 4.13g NaF per lOOg water) and KF (15.4M, or 89.8g KF 
per lOOg water) in water at 20°C. In comparison, the saturation solubilities for the chloride 



and bromide salts of these three cations are not so different. [36| These facts suggest that 
such differences are most possibly due to the pairing or association of the related anion and 
cation, since no single ion effect can explain their differences. In the Table Ull the k value 
of the NaF clearly stood out of the other five salts including KF. Similar to the explanation 
for their bulk solubilities, this is also most likely due to the Na^ and F~ ion pairing or 
association, instead of the property of the F~ anion alone. 

Ion pairing is most likely to be found in aqueous solutions at high bulk concentration. For 
example, a 5M monovalent salt solution has less than six water molecules for each ion in the 



solution. In order to fully solvate an ion in the solvation shell, usually six water molecules 
72, l75| At such high bulk concentration as 5M the cations and anions are 



are needed. 
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very close to each other. However, our SHG data of the five salts KF, KCl, KBr, NaBr and 
NaCl did not show the same ion pairing or association effects as the NaF salt. Therefore, 
ion pairing or ion association can not be simply defined by the average distance between the 
cations and anions even at high bulk concentrations. 

In the end, we have to realize that all the effects as measured with the non- resonant SHG, 
i.e. the increase of the average thickness of the interfacial water layer, as well as the shift 
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of the averaged orientational angle of the surface water molecules, are really small. The 
increase of the thickness of the surface water layer is most significant for the NaBr and KBr 
aqueous solutions. The increase is only less than 30% for the 5M NaBr salt solution (Figure 
[3]). This increase is much smaller than the more than 100% thickness increase value for the 
2.0M (0.036 molar fraction) NaBr solution surface estimated using the relative Raman, IR 



and SFG intensities in the literature. 



27( 1 Increase of the surface water layer thickness requires 



more ions to be solvated in the surface water layer. Otherwise the surface concentration of 
the ions could not keep up with the bulk ion concentration. Therefore, in addition to the 
quantitative treatment in the non-resonant SHG, we surmise that the slight increases of the 
average water surface layer thickness as reported here are more acceptable. Furthermore, 
the average orientational parameter D for the surface water molecules also changed slightly 
for all the salt solutions studied (Figure [2]). All these suggested that the ion effects on the 
surface water molecules were subtle effects. However, no matter how weak these effects 
were, the presence and direct interaction of the ions to the surface water molecules have 
been firmly established. These conclusions may add vindications to the recent interests on 
the specific ion effects at the air /water surfaces. 

In the experimental studies on the reactive and non-reactive molecular beam scattering 
from the salty liquid solution surfaces, such as from the glycerol salt solution surfaces, various 



cation and anion effects have been observed. 
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8l| Since molecular beam scattering 



would be difficult to be implemented for the more volatile liquids than the glycerol, such 
as water, it would be worthwhile to investigate the ion solvation and interactions at these 
non-aqueous liquid surfaces with the nonlinear optical techniques and to compare the results. 

In summary, the specific ion effects on the surface water molecules at the electrolyte 
aqueous solution surfaces can be more complex than expected. Quantitative surface specific 
measurement such as the non-resonant SHG can be employed to help establishment of the 
detailed knowledge about these surfaces. Since the observed specific cation and anion effects 
on the surface water molecules are in the comparable magnitude, theoretical modeling and 
simulation efforts are called upon to understand these unexpected cation effects for the Na"*" 
and K+ ions, as well as the unexpected anion effects for the F~ ions, which are usually 
considered to be repelled from the aqueous solution surface. 
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V. CONCLUSION 



In this report, using the procedures for sample control and quantitative SHG data anal- 
ysis, we studied the non-resonant SHG from the water molecules at the aqueous solution 
surfaces of the following monovalent salts: NaF, NaCl, NaBr, KF, KCl, and KBr. In addi- 
tion to the F^, Cl^ and Br~ anion effects as reported in our previous report, the specific Na^ 
and cation effects on the interfacial water molecules are also clearly and quantitatively 
identified for the first time. 

Through quantitative polarization analysis of the non-resonant SHG data, we found that 
the orientational order parameter D of the water molecules at the interfacial region changed 
in different direction for the NaX (X=F, CI, or Br) series versus the KX series salts. The 
average orientation angle of the dipole vector of the water molecules at the NaX aqueous 
solution surfaces changed slightly to the smaller tilt angle from the surface normal in compar- 
ison to that of the neat air/water interface; while for the KX series, the average orientation 
angle changed slightly to the larger tilt angle. This seems to suggest that the Na"*" cation is 
more "structure-making", while the K"*" is more "structure-breaking". The relative ability 
for the different salts to increase the average thickness of the interface water layer was also 
quantitatively determined for the six salts. We found that these values can be factorized 
into simple additive values for the NaCl, NaBr, KF, KCl, and KBr salts except for the NaF. 
This indicated that ion pairing or association effects for the Na"*" and F~ ions should be 
much stronger than the cation and anion in the other five salts. 

Specific Na"*" and K"*" cation effects can also be observed in the SFG-VS spectra of the 



interfacial water molecules at the NaF and KF aqueous solution surfaces. [62] These results on 
the specific cation and anion effects suggest that the electrolyte aqueous solution surfaces can 
be more complex than the currently prevalent theoretical and experimental understandings. 
Since the observed specific cation and anion effects on the surface water molecules are 
in comparable order of magnitudes, theoretical modeling and simulation efforts are called 
upon to understand these unexpected cation effects for the Na"*" and K"*" ions, as well as the 
unexpected anion effects for the F~ ions, which are usually considered to be repelled from 
the aqueous solution surface. 
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Table 1 The results for the XzxxIXxzx-, Xzzz/Xxzx, R-, D, 6 (assuming 5 distribution) and 
for the neat water, KF, KCl and KBr salt solution interfaces. 

Table 2 The relative ability for each salt to increase the thickness of the interfacial water 
is determined by fitting the slopes in the Figure [31 The six salts followed the 
following order: KBr > NaBr > KCl > NaCl ~ NaF > KF. 
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TABLE I: The results for the XzxxjXxzx^ Xzzz/Xxzx-, R, D, 6 (assuming 5 distribution) and 
the neat water, KF, KCl and KBr salt solution interfaces. 



KX Concentration/M Xzxx/Xxzx Xzzz/Xxzx R B 9 Ng 


H2O 


0.0 


0.29±0.03 2.53±0.06 0.69±0.02 1.71±0.02 40.1±0.4 1 


KF 


0.5 


30±0 08 2 51 ±0 06 BQiO 02 1 72±0 02 40 8±0 4 1 OliO 02 




1.0 


30±0 03 2 57±0 06 69±0 02 1 72±0 02 40 3ib0 4 1 01±0 02 




2.0 


31ib0 03 2 58±0 06 70±0 02 1 72±0 02 40 3±0 4 1 02±0 02 




3.0 


31±0 03 2 58±0 06 70±0 02 1 72±0 02 40 3±0 4 1 04±0 02 




4.0 


32±0 03 2 62±0 06 70±0 02 1 72±0 02 40 3±0 4 1 05±0 02 




5.0 


0.33±0.03 2.60±0.06 0.71±0.02 1.72±0.02 40.3±0.4 1.07±0.02 


KCl 


0.5 


29±0 03 2 52ib0 06 69±0 02 1 71±0 02 40 1±0 4 1 02±0 02 




1.0 


0.29±0.03 2.57±0.06 0.69±0.02 1.72±0.02 40.3±0.4 1.03±0.02 




2.0 


0.28±0.03 2.53±0.06 0.68±0.02 1.71±0.02 40.1±0.4 1.08±0.02 




2.5 


0.28±0.03 2.52±0.06 0.68±0.02 1.71±0.02 40.1±0.4 1.10±0.02 




3.0 


0.29±0.03 2.56±0.06 0.69±0.02 1.72±0.02 40.3±0.4 1.12±0.02 




4.0 


0.28±0.03 2.54±0.06 0.68±0.02 1.72±0.02 40.3±0.4 1.16±0.02 


KBr 


0.5 


0.29±0.03 2.54±0.06 0.69±0.02 1.71±0.02 40.1±0.4 1.05±0.02 




1.0 


0.28±0.03 2.54±0.06 0.68±0.02 1.72±0.02 40.3±0.4 1.08±0.02 




2.0 


0.27±0.03 2.59±0.06 0.68±0.02 1.73±0.02 40.5±0.4 1.17±0.02 




3.0 


0.27±0.03 2.61±0.06 0.68±0.02 1.73±0.02 40.5±0.4 1.25±0.02 




4.0 


0.27±0.03 2.57±0.06 0.68±0.02 1.72±0.02 40.3±0.4 1.31±0.02 
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TABLE II: The relative ability for each salt to increase the thickness of the interfacial water is 
determined by fitting the slopes in the Figure [3l The six salts followed the following order: KBr 
> NaBr > KCl > NaCl ~ NaF > KF. 

Salts KBr NaBr KCl NaCl NaF KF 

Slope (M-i) O.OSOibO.OOl 0.069±0.001 0.040±0.001 0.028±0.001 0.029±0.002 0.013±0.001 
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Fig. 1 Polarization and concentration dependent SHG (400nm) data: a. left panel, KF, 
KCl and KBr solution surfaces; b. right panel, NaF, NaCl and NaBr solution 
surfaces. These data were obtained by fitting the s and p detection data measured 
for each concentration and each salt solution. The procedure used was fully 
described previously, Q, Q and it can ensure the accuracy of the data for each 
polarization combination. The 45°-in/s-out data need to be scaled up by the 
factor of 1.23 because of the different detection efficiency for the s and p polarized 
SHG signal at 400nm in the detection setup in this experiment. The data for 
NaF, NaCl and NaBr are identical to the data in our previous article, [l^ and the 
scale is changed to be the same as the KF, KCl and KBr data for better visual 
comparison. 

Fig. 2 The orientational parameters D (top) and the apparent orientational angles 9 
(bottom) for the KF (solid triangle), KCl(solid circle), KBr (solid square), NaF 
(open triangle), NaCl(open circle) and NaBr (open square) solution interfaces. 
The data for KF, KCl and KBr are from Table [H The data for NaF, NaCl and 
NaBr were from the Table 1 in the previously published paper. It is clear that 
the D values for the neat water and low bulk concentration for the sodium halides 
and the potassium halides data series were well within the experimental error; 
while at higher concentrations the D values for the two series went to different 
directions. 

Fig. 3 The relative interfacial density of the water molecules at different salt bulk con- 
centration: KF (solid triangle), KCl(solid circle), KBr (solid square), NaF (open 
triangle), NaCl(open circle) and NaBr (open square). The data for KF, KCl and 
KBr are from Table [B The data for NaF, NaCl and NaBr were the same as in 
the previously published work. 1^ 

Fig. 4 The polarization dependent per molecule hyperpolarizabilities (Xpp.o, X45°s,05 and 
Xspfl) are simulated with Eqs J2.2l and Eqs l2.4[ In calculating the Fresnel fac- 
tors, the effective surface dielectric constant was calculated according Zhuang 



et al.^s model. [38| The experimental configuration is in the refiection geometry 
with P=70°. The values of R={Pcaa + Pcbb) / {Paca + Pbcb)=0 .69 and {Paca + Pbcb)=f 
are used. In the simulation, the angular distribution function was assumed a 
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Gaussian function f{9 — 9o) — :^/^e ^0)^/20-2 rpj^^ values for 5 distribution 
function (cr = 0°) and for FWHM=30° ((7=12.7°) are both plotted. The vertical 
dash line is at the average orientational angle (6) of the neat air /water interface 
which was determined at 40.1°. The horizontal dash line is drawn at Xe//,o=0. 

Fig. 5 Illustration of the two types of interfacial water structure as well as the specific 
cation effects on the average orientation of the non-straddle type surface water 
molecules. There may exist different type of configurations for the non-straddle 
type water molecules. Here only their averaged orientation of these hydrogen 
bonded water molecules is representatively illustrated. 
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FIG. 1: Polarization and concentration dependent SHG (400nm) data: a. left panel, KF, KCl 
and KBr solution surfaces; b. right panel, NaF, NaCl and NaBr solution surfaces. These data 
were obtained by fitting the s and p detection data measured for each concentration and each salt 
solution. The procedure used was fully described previously, and it can ensure the accuracy 

of the data for each polarization combination. The 45°-in/s-out data need to be scaled up by the 
factor of 1.23 because of the different detection efficiency for the s and p polarized SHG signal at 
400nm in the detection setup in this experiment. The data for NaF, NaCl and NaBr are identical 
to the data in our previous article, [14] and the scale is changed to be the same as the KF, KCl and 
KBr data for better visual comparison. 
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FIG. 2: The orientational parameters D (top) and the apparent orientational angles 6 (bottom) 
for the KF (sohd triangle), KCl(solid circle), KBr (solid square), NaF (open triangle), NaCl(open 
circle) and NaBr (open square) solution interfaces. The data for KF, KCl and KBr are from Table 
m The data for NaF, NaCl and NaBr were from the Table 1 in the previously published paper, [l^ 
It is clear that the D values for the neat water and low bulk concentration for the sodium halides 
and the potassium halides data series were well within the experimental error; while at higher 
concentrations the D values for the two series went to different directions. 
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FIG. 3: The relative interfacial density of the water molecules at different salt bulk concentration: 
KF (solid triangle), KCl(solid circle), KBr (solid square), NaF (open triangle), NaCl(open circle) 
and NaBr (open square). The data for KF, KCl and KBr are from Table [H The data for NaF, 



NaCl and NaBr were the same as in the previously published work. 
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FIG. 4: The polarization dependent per molecule hyperpolarizabilities (Xpp,0) Xa5°s,o, and Xsp,o) are 
simulated with Eqs l2.2l and Eqs l2.4[ In calculating the Presnel factors, the effective surface dielectric 
constant was calculated according Zhuang et aZ.'s model. [Q] The experimental configuration is in 
the reflection geometry with 13=70°. The values of R={f3caa + Pcbb) / {f^aca + Pbcb)=0.69 and {(3aca + 
Pbcb)=^ are used. In the simulation, the angular distribution function was assumed a Gaussian 
function f{d — 9q) = g-(6'-6'o)^/2o-^ ^ The values for 6 distribution function (cr = 0°) and for 
FWHM=30° ((7=12.7°) are both plotted. The vertical dash line is at the average orientational 
angle {0) of the neat air/water interface which was determined at 40.1°. The horizontal dash line 
is drawn at Xe//,o=0. 



32 



Air 



Interfacial region 



NaX, (X=F, CI or Br) ^ I 
Slightly less tilted 
Slightly thicl<er surface 



Bulii 



30-35° 




Free 0-H ( 3700cm-i ) 
(Straddle type HjO) 
SFG detectable 

Hydrogen bonded HjO 
(Non-straddle type HjO) 
SHG & SFG detectable 



f Slightly more tilted 

Slightly thicker surface 



FIG. 5: Illustration of the two types of interfacial water structure as well as the specific cation 
effects on the average orientation of the non-straddle type surface water molecules. There may 
exist different type of configurations for the non-straddle type water molecules. Here only their 
averaged orientation of these hydrogen bonded water molecules is representatively illustrated. 
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